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Besides the growing interest in old concepts like temperature and entropy at the 
nanoscale, theories of relaxation and transport have recently regained a lot of attention. 
With the electronic circuits and computer chips getting smaller and smaller, a fresh look 
should be appropriate on the equilibrium and nonequilibrium thermodynamics at small 
length scales far below the thermodynamic limit, i.e. on the theoretical understanding of 
original macroscopic processes, e.g. transport of energy, heat, charge, mass, magnetiza- 
tion etc. Only from the foundations of a theory its limits of applicability may be inferred. 
This review tries to give an overview of the background and recent developments in the 
field of nonequilibrium quantum thermodynamics focusing on the transport of heat in 
small quantum systems. 
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1. Introduction 

How is energy or heat transported through a material system on a microscopic level 
of description? Despite its long history this old question has remained exciting even 
today, since a satisfactory answer has not yet been found. The large number of 
different publications on the topic of heat transport indicates both an increasing 
interest in the problem and its nontrivial character, emphasized by the title of a 
recent publication by Bonetto^ et al. "Fourier's Law: A Challenge to Theorists". 
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The review at hand tries to give an introduction to the phenomenon of transport 
of heat from the viewpoint of various theoretical approaches, before turning to 
some new perspectives on this old topic. In such nonequilibrium phenomena as 
heat conduction, the intensive parameters of different parts of a system are no 
longer equal, i.e., there is no global equilibrium established inside the whole system. 
However, despite of this fact, in a macroscopically small but microscopically large- 
part of the system, at least near the global equilibrium, one finds an equilibrium 
state slightly different from region to region. One may think, for example, of a 
material bar coupled at both ends to heat baths of different temperature. As noted 
by Fourier— 

When heat is unequally distributed among the different parts of a solid 
mass, it tends to attain equilibrium, and passes slowly from the parts which 
are more heated to those which are less; and at the same time it is dissipated 
at the surface, and lost in the medium or in the void. 

In this sort of experiment one typically finds a constant temperature gradient within 
the material. Thus, there is no global temperature defined in the whole system, but 
in a small enough part one could approximately describe the system by a canonical 
equilibrium state with local temperature. 

The first detailed theoretical as well as experimental investigation of heat trans- 
port through solids was carried out by Joseph Fourier'^ about 1807. The subsequent 
Paragraph is dedicated to him introducing the famous Fourier's Law. In the follow- 
ing we will review some of the main approaches to the famous law on microscopic 
level as well as their respective shortcomings, before we give an overview of recent 
developments in the heat conduction research. The second Section will introduce 
some main ingredients of nonequilibrium physics followed by the analysis of heat 
conduction in small quantum systems by the quantum master equation. The last 
part of this review presents a recently developed theory, the so-called Hilbert Space 
Average Method (HAM), to study the relaxation into equilibrium. This method 
allows to approach thermal transport from a completely different viewpoint. 

2. Overview 

2.1. Fourier's Law 

Almost two hundred years ago Fourier conjectured that temperature (or as we know 
today: energy) tends to diffuse^ through solids once close enough to equilibrium. 
In his work he considered the transport of heat on a phenomenological respectively 
macroscopic level, defining a partial differential equation for the heat transport 
through some material and trying to solve it by mathematical methods, e.g., the 
Fourier series expansion. Nowadays, Fourier is mainly known for his marvelous 
mathematical research, especially the celebrated Fourier analysis. 

Considering a situation with a position and time dependent temperature profile 
T{r,t) within some material, Fourier describes the problem by an equation for the 
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internal energy density u{T{r, t)). The change of this quantity u{T{r, t)) should be 
proportional to the gradient of the thermodynamic force - here the temperature 
gradient 'VT{r,t), 

^w(r(r,i)) =«V- VT(r,i) , (1) 

with K being the heat conductivity. As the energy density is a function of temper- 
ature one can also state that 

and by using the definition of the specific heat capacity 

c(T) = ^ (3) 

one finds the heat conduction equation 

|T(r,.) = -^AT(.,.). (4) 

This diffusive behavior is slightly reformulated by applying the continuity equation 
for the energy density 

l + V.J.O, (5) 

where J represents the heat current in the system. Plugging ([5]) into one finds 

V V (KVT(r,i)) . (6) 

From these considerations one immediately gets what Fourier's law states close to 
thermal equilibrium: a proportionality between the heat current and the tempera- 
ture gradient 

J = -KVT{r,t). (7) 

Furthermore by thinking of the temperature T being a function of the internal 
energy density u we could state for the energy current 

dT K 

J = -nVTiuir, t)) = -n—Vu = Vu(r, t) , (8) 

' ou c 

where J is now the energy current density. 

Despite the ubiquitous occurrence of the phenomenon of heat conduction in our 
everyday experience, its explanation on the basis of some reversible microscopic 
dynamics remains a serious problem. The quantity of interest here is the material 
property k - the heat conductivity. In the last two centuries there have been several 
different microscopic approaches to Fourier's Law from completely different back- 
grounds. The goal of all of those efforts was to give a thorough and comprehensive 
picture of heat transport from the viewpoint of a basic microscopic theory. 
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2.2. Kinetic Gas Theory for Phonons 

One of the first attempts to find a microscopic foundation for the heat conductivity 
of an insulating sohd was due to Debye'^'^. He considered the phonons inside such a 
material as classical particles using results of the kinetic gas theory to describe the 
heat transport based on phonons. 

The particle current in one direction of a crystal should be proportional to the 
density rj of the particles and their mean velocity v. A single phonon considered 
as a classical particle moving from a region with temperature T + AT to a region 
with temperature T looses the energy CAT with C being the heat capacity of the 
respective particle. With the internal temperature gradient VT along the free path 
I of the phonon we find the temperature difference to be 

AT = Ivr. (9) 

In such a situation the energy current is just proportional to the particle current 
times the transported energy and thus 

J (XTjvCAT ocTjvClVT . (10) 

Reformulating this result by the heat capacity per volume c = rjC and comparing 
it to Fourier's Law ^ we find for the heat conductivity coefficient 

K oc cvl . (11) 

As the free path of the phonons enters the above microscopic formula an infinite 
heat conductivity would result in case of no scattering in the solid. Why should 
phonons have a finite free path inside a material? Debye"* and others have mainly 
taken two processes into account which limit the free path of the phonons, the 
scattering at impurities (primarily in amorphous materials) and the phonon-phonon 
interaction in the solid (perfect crystal, lacking any impurities). However, if the 
lattice was harmonic there would be no scattering between phonons and therefore 
no finite heat conductivity. 

To include anharmonicity in the coupling of the atoms in the solid is a serious 
theoretical problem. Debye* proposed that the Brownian motion of the atoms might 
produce an irregularity in the lattice and thus give rise to phonon scattering in an 
otherwise perfect crystal. On this basis he was able to account for the mean free path 
of the phonons, finding it anti-proportional to the temperature of the solid. In this 
way he was able to derive an approximate formula for a finite heat conductivity— 
from a microscopic point of view. 

2.3. Peierls-Boltzmann Equation 

A more detailed approach to the question of a microscopic foundation of Fourier's 
Law was introduced by Peierls- ' - in 1929 applying the Boltzmann equation to the 
phonons in the solid. Peierls mainly considered the mentioned anharmonicities in 
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the atom coupling, which are the dominant effect underlying phonon scattering in 
pure crystals of not too small size and at not too low temperatures. 

Treating the anharmonic part of the interaction as a perturbation, Peierls was 
able to find conditions from Fermi's Golden Rule for the possible scattering pro- 
cesses: the energy conservation and the momentum conservation 



where (cji, fci) and (^2,^2) are the two incoming phonons, (103, k^) is the outgoing 
one, and g is the appropriate reciprocal lattice vector—. All processes with g = are 
not suitable for a thermalization of the phonon gas, since then the whole momentum 
of all phonons is conserved, the conductivity would again be infinite (all these 
processes are called normal or N-processes). Of vital importance for the existence 
of a thermal resistance are the scattering processes with g ^ called Umklapprozess. 

Considering now very low temperatures {T < Q, Q Debye temperature) where 
most of the two incoming thermally activated phonons have not enough energy to 
fulfill the above scattering conditions, there is again no possible umklapprocess. 
According to the Boltzmann distribution of phonon energies we expect that for 
low temperatures the number of high-energy phonons decreases by an exponential 
g-e/T^ Thus, the conductivity must be proportional to that exponential factor in 
the low temperature regime. 

Since phonons are also scattered at impurities or boundaries of the crystal, one 
has to consider more scattering processes in order to get a complete picture of heat 
conduction in solid&^>^'^ '^''i^^ . These additional processes are very important for low 
temperatures, especially, since in this region the umklapprocesses die out rapidly. 

To describe the phonons inside the crystal Peierls essentially proposed a Boltz- 
mann equation replacing classical gas particles by quantized quasi particles - the 
phonons. Furthermore he assumed classical transition probabilities between differ- 
ent phonon modes as obtained from Fermi's Golden Rule. His theory succeeds in 
explaining the basic transport properties of solids and allows to derive statements 
about the conductivity of insulating crystals. The equation itself is very hard to 
solve, but, nevertheless, one can obtain some tendencies in the temperature depen- 
dence of the conductivity. Peierls^ states: 

The only prediction that would seem to follow with certainty is that the 
law must be intermediate between T^^ and T^^. 

Here the different dependencies between 1/T and follow from third order 

respectively forth order anharmonicitics in the interaction. 

However, the Peierls-Boltzmann theory faces conceptual shortcomings: quan- 
tized normal modes of a many particle system are described classically, i.e., they 
are well localized in configuration as well as momentum space and show no disper- 
sion—. This becomes a serious problem at least in the limit of systems consisting of 
only a few particles at all. Already Debye's approach to heat conduction deals with 



UJl + UJ2 



fci + A;2 = fcs + , 71 G Z , 



^3 , 



(12) 
(13) 
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a rather classical idea of the quasi-particles - thinking of bouncing balls like gas par- 
ticles in a box instead of quantum mechanical objects. Furthermore, in order to ex- 
ploit Fermi's Golden Rule, the actual quantum state of a phonon mode is discarded 
keeping only the mean occupation number. Due to the neglect of any phases, this is 
called the random phase approximation and replaces Boltzmann's Stosszahlansatz 
(see Peierls^). Strictly speaking these both lack justification, though, and there are 
no clear criteria for the applicability of the theory. Additionally, the application 
of the theoretical concept to a concrete coupling model is very complicated and 
Peierlsi himself states: 

This seems a poor return for a long discussion, but progress beyond this 
stage is difficult, unless we could construct a dispersion law which was 
simple enough to allow us to list the solutions for the possible phonon 
collisions explicitly, and yet realistic enough to give the right kind of collision 
including those of the Umklapp type. 

Even today there are doubts about the importance of umklapprocesse s ^^'^'^ at all 
for a finite conductivity in solid states. 

2.4. Linear Response Theory 

Another powerful technique within the field of heat transport is the Green-Kubo 
formula. Derived on the basis of linear response theory it has originally been for- 
mulated for electrical transport i^^iiSiie ^ Here the external force, the electric field, is 
considered as the perturbation of the system. In this context the current is viewed 
as the response to this external perturbative electrical potential which can be ex- 
pressed as a part of the Hamiltonian of the system—. 

2.4.1. Currents, Forces and Onsager Relation 

In a linear response theory, currents Ji are linear function of external forces Xj 
defining the matrix of response coefficients l-ij as 



From the Onsager relation we know that the coefficient matrix has to be symmetrical 



The choice for the current operators and the corresponding forces is ambiguous. 
But it turns out that the Onsager relation is not valid for all choices. A special 
one for which the Onsager relation is valid was proposed by de Groot^^, who chose 
currents and forces implicitly defined by the expression for the increase of entropy 



since a nonequilibrium process should be associated with a net entropy production. 




(14) 



■ 15 .17.18 




(15) 
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Let us consider a situation with two different fluxes, a particle current Ji due to 
an external force Xi = — 7^V/i with the generalized potential /i (contains chemical 
as well as electrical potentials) and a heat current J2 due to the force X2 = Vy. 
Then ([14]) reads 

Ji = -^ZiiVm + Zi2V^ , (16) 

J2 = -^Z2lV/i + Z22V^ . (17) 

For a vanishing particle current Ji = it follows from (jl6p 

iZnVM-Zi2vi (18) 
and, using Onsager's relation and reformulating pT)) with ([18]), for the heat current 

J2=(Z22-^)V1. (19) 

Putting aside any further perturbations like chemical or electrical potentials (/i = 0) 
we obtain 

J2 = Z22vi - Z22 (-^V^) ■ (20) 
Comparing this result to Fourier's Law ([7]) one finds for the heat conductivity 

«=^- (21) 

For a microscopic theory of the heat conductivity it remains to compute the response 
coefficients Z22- 



2.4.2. Transport Coefficient 

In case of an electric perturbation of a system Kubo was able to derive a formula 
for the conductivity tensor cry (uj) by first order perturbation theory. In so doing 
he considered the external electric potential given as an operator part within the 
Hamiltonian of the system. The result of this consideration is basically a current- 
current autocorrelation function— 

a^j{uj)= dte'"* / dXTi {Jj{-t - ihX) po} , (22) 
Jo Jo 

with the frequency to of the perturbation and po being the equilibrium state of the 
system at inverse temperature /3. The current operators are here in the Heisenberg 
picture, i.e., for a system with an unperturbed Hamiltonian H they are defined as 

J(x) := e*^^ Je"-^^^ . (23) 

Since the above Kubo formula essentially consists of a current-current auto- 
correlation, it may ad hoc be transfered to heat transport simply by replacing the 
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electrical current by the heat current. However, the justification of this replacement 
remains a conceptual problem since there is no way of expressing a temperature 
gradient in terms of an addend to the Hamiltonian which has been no problem for the 
electric case. Remarkably enough, Kubo himself commented on that replacement 
in a rather critical way: 

Therefore, the treatment developed ...[above] does not to be directly appli- 
cable to nonequilibrium states produced by such thermal forces, in order 
to obtain explicit formulas expressing the responses to thermal forces. It is 
generally accepted, however, that such formulas exist and are of the same 
form as those for responses to mechanical disturbances. 

Nevertheless, Luttinger^^ proposed in 1964 a formula for the response coefficient 
for thermal perturbations on the basis of Kubo's theory 

Z22=/ dte-'* dXTr{J{-t-ihX)Jpo} . (24) 
Jo Jo 

Additionally, one needs the limit of a; ^ here since a frequency dependent thermal 
perturbation seems even more suspicious than the questionable analogy to electric 
perturbations. The exponential function here with s > is a switch-on function of 
the perturbation from the infinite past, to guarantee the convergence of the integral. 

There are some further discussions about the divergence of the response coeffi- 
cient for w ^ and the respective weight of the delta-function (Drude- weight) as 
well as the regular part inside the above formula. But these discussions are definitely 
beyond the scope of this text and therefore we would like to refer to the literature 
(see e.g. the work by Heidrich-Meisner— ) . 

Since there is no detailed derivation based on assumptions or preconditions, it 
is hard to state the limits of applicability of the formula. Despite all insufficiencies 
of this approach, it has become a widely employed techniqu e and it allows 
for a straightforward application to any system, once it is partially diagonalized. 

2.5. Survey of Recent Developments 

During the past decade the established field of "heat transport" has regained at- 
tention in the physics community. This is not only a consequence of the conceptual 
difficulties of the above described standard microscopic approaches - the Peierls- 
Boltzmann theory and the Green-Kubo formula - but is also a result of some new 
and, in part, technical developments. 

In a time where electronic circuits and computers are getting smaller and smaller, 
the knowledge about thermodynamical properties on the nanoscale becomes more 
and more important. Especially, the quantum limits to some standard thermody- 
namical concepts may eventually lead to a revolution not only in the theoretical 
and experimental understanding of the dynamics of small systems, but also in their 
technical applications. Thus, there is an increasing interest in experimental as well 
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as theoretical advancements concerning heat transport in nanoscopic devices. All 
together this has led to a revival of old, but to a certain degree unsolved ques- 
tions and to a new vital discussion on heat conduction, especially in small quantum 
systems. 

Unfortunately, the microscopic foundation of Fourier's Law seems to be very 
complicated and far from being trivial as stated by the above mentioned article 
by Bonetto^'^ . The numerous difficulties tempted the authors to promise a bottle of 
wine for a proper microscopic theory. Another statement about the actual status of 
heat conduction research appeared recently in an article of Buchanan^"* in Nature 
Physics: "No one has yet managed to derive Fourier's Law truly from fundamental 
principles" . However, we hope to convince the reader of this text that there are now 
several ideas of how to approach normal heat conduction and thus Fourier's Law 
from quantum mechanics directly. 

The following selection of heat conduction research can hardly do justice to 
respective approaches nor claims completeness. It is only meant to show the great 
variety of different ideas and their links to our work. 

2.5.1. Heat Conduction in Classical Systems 

Since the work at hand is mainly based on quantum mechanical ideas, we only 
briefly comment on classical models of heat conduction. We refer the interested 
reader to the excellent review article by Lepri— et al. , which summarizes the central 
techniques and the main results of heat conduction in classical mechanics. 

In the classical domain it seems to be largely accepted that normal transport 
requires chaotic microscopic dynamics^^'^^ whereas ballistic transport is typically to 
be found in completely integrable systems following a regular dynamics. Therefore, 
one finds, e.g., a normal transport in a Lorentz gas^L^ model which is strongly 
chaotic (of course, the heat transport is strictly connected to particle diffusion, here) . 
However, there have also been successful attempts to observe normal transport in 
the absence of exponential instability—, i.e. without chaos. 

In view of the many different but often strange and unrealistic model systems 
under consideration one gets the impression that it is really hard to find the desired 
diffusive behavior. This contradicts our all day experience, since normal diffusion 
clearly is the typical property rather than ballistic behavior. Nevertheless, people 
have proposed several "spring and ball" models combined with classical free particle 
models and different masses etc. ( "ding-a-ling" -model— , "ding-a-dong" -model—) as 
well as chains of nonlinear oscillator s siiS^^sa, 34^35, 36 finally obtain normal trans- 
port behavior. 

2.5.2. Experimental Investigations 

Recently, there have been some efforts to measure the heat conductivity of a single 
carbon nanotube. Besides some evidences of ballistic phonon transport in carbon 
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nanotubes^i^, there are measurements on the length and temperature dependence 
of the heat conductivity ^9 ^ too. 

Furthermore, there are several investigations on the thermal conductance of 
magnetic systems. Those materials are insulators, for which the heat transport is not 
dominated by electrons. Inside such materials one frequently finds chains, ladders or 
quasi two dimensional spin structures. Thus, besides the standard phonon induced 
heat transport in insulators, also magnons could be involved. A characteristic of 
such a behavior is the large anisotropy between the thermal conductivity measured 
parallel to spin structures and perpendicular to them. Furthermore, the absolute 
value of the heat conductivity in some insulating compounds is extremely high, and 
comparable with the conductivity of some metals. This has led to the assumption 
of ballistic magnon transport in such materials. Recent measurements have been 
carried out by Sologubenko^MLii^d^ et al. and Hess^i^i^ et al. For a nice overview 
on experimental results in different materials cf. the work of Heidrich-Meisner— . 

2.5.3. Kuho-Formula 

To support the experimental results mentioned in the last Section there are several 
attempts to evaluate the thermal Kubo-formula for spin structures. Inspite of the 
above discussed insufficiencies of the Kubo-formula, it is, nevertheless, a frequently 
used technique to account for the thermal conductivity of a system. Furthermore, it 
allows for a discussion of the emergence of regular transport, scaling properties of the 
conductivity and temperature dependencies. Besides some analytical attempts for 
integrable system&2ii^, there are numerous numerical one s 20,22,23,48,49,50,51 ^ com- 
puting the Kubo-formula for several concrete spin systems. The main quantity is the 
Drude weight, i.e. the conductivity for frequency going to zero. A divergent weight 
indicates ballistic transport. One promising method evaluates the whole frequency 
dependent conductivity first, a set of delta-peaks with different weight for finite 
systems. Sorting those peaks into bins and estimating the value for frequency to 
zero, results in the desired conductivity—!^. 

The big advantage of the Kubo formula is its computability for a concrete model 
after having partially diagonalized its Hamiltonian. However, it comes without any 
criteria for its applicability, thus, potentially leading to wrong results. According to 
a new theory— it seems to be possible to find the celebrated formula for thermal 
transport by first principles. 

2.5.4. Reservoir Coupling 

Essentially, the Kubo formula is based on a linear response theory. But, in order 
to derive the proper equation, the external perturbation has to be formulated as 
a potential part in the Hamiltonian of the system. As discussed above this is not 
possible for thermal perturbations. To overcome this problem, such Kubo-scenarios 
have recently been transfered from Hilbert- to Liouville space, where temperature 
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gradients may indeed be formulated in terms of operators^*. In this approach several 
heat reservoirs are coupled to the system described by a quantum master equation. 
Numerically, the method reveals normal heat transport ss^se, 57^58, 59^60 final 
stationary local equilibrium state, already in surprisingly small quantum systems. 
Furthermore, there are some attempts to solve such bath scenarios analytically, e.g. 
for the Ising model—. 

Another approach based on the coupling the system to different heat baths leads 
to the celebrated Landauer-Biittiker formula Here, the system is considered as a 
junction, e.g. a small molecule, which transmits or reflects incoming modes from the 
two reservoirs of different temperature. The heat flow is thus defined only by the 
transmitted modes. Again the model reveals both normal and ballistic transport in 
dependence of the type of the junction— iHi^. However, since the system does not 
contain an internal geometrical substructure, there is no room for getting Fourier's 
Law inside the system proper. 

Instead of the coupling to reservoirs of different temperature in Liouville space, 
there are proposals to impose a current inside the system by adding the standard 
current operator— i^i^ to the Hamiltonian of the system. Finally, one observes 
again the quasi stationary state in the system finding domains of normal respectively 
ballistic transport in spin chains. 

Finally, there are efforts to study the heat conducting behavior of large spin 
systems by the new numerical method of time dependent density matrix renormal- 
ization groupS-^ (tDMRG). 

2.5.5. Quantum Chaos 

Just as in the classical domain there are some investigations on the correlation be- 
tween regular transport and the onset of quantum chaos, too. Since the Schrodinger 
equation is a linear equation of motion, the chaos debate in quantum mechanics and 
the onset of "quantum chaos" itself is an open question up to now. However, Mejia- 
Monasterio^° et al. have recently found some evidences of chaotic behavior of a spin 
chain model under a heat conducting bath coupling. By passing from ballistic to 
regular transport as a single system parameter is changed, the system switches from 
an integrable model into a chaotic one. This onset of quantum chaos is observed 
by the change in the level distribution form a Poissonian to a Wigncr-Dyson type 
distribution^^. There have been attempts to find correlations between the change 
of the level statistic of a system and its transport behavior''^. 

3. Description of Nonequilibrium Scenarios 
3.1. Topological Structure 

Quantum mechanical systems are described by a Hamilton operator H defined in 
the respective Hilbert space Ti of the system. The class of systems considered here 
consists of those made up of several subunits, n-levels each, coupled by different 
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types of interactions. Therefore the Hamiltonian of the complete system is described 
by a local part for each subimit /i, Hiodp-), and an interaction, Hint- In the following 
we concentrate on one-dimensional systems (chains) of N subunits and on a next 
neighbor interaction, only. Accordingly, the complete n^-dimensional Hamiltonian 
reads 

N N-1 

i?-^Aoc(/i)+A^7?int(Ai,M+l), (25) 
/^— 1 M— 1 

with A denoting the coupling strength. This coupling should be weak in order to 
guarantee the survival of local features. 

A simple example for this model would be a spin chain. In terms of Pauli oper- 
ators (Ti (i = x,y,z) the local Hamiltonian of a subunit /i with an energy splitting 
AE can then be written as 

Aoc(/i) = — <7.(m) (26) 
while the next neighbor interaction reads 

Ant(M,M+l)= X! Ci(Ti(^)(»(7i(^ -I- 1) . (27) 

This is the Heisenberg interaction for Cx — Cy — Cz — const, and the so-called 
Forster coupling for = 0. 

3.2. The Heat Current 

For a system consisting of several subunits and a Hamiltonian like ([25|) . we expect 
most of the energy to be concentrated in the local part of the Hamiltonian of the 
system, only a small amount being in the interaction (weak coupling limit). To get 
an operator for the current^i2ii^i2^i2i between two adjacent subunits in the system, 
we consider the time evolution of the local energy operator given by the Heisenberg 
equation of motion 

d - , , i ~ - , dH]„c(u) 

-i7ioc(A^) = j: [H, H^M] + Q^' ■ (28) 

Since here and in the following the local Hamiltonian does not explicitly depend 
on time, the last term vanishes. Plugging in the complete Hamiltonian (j25p and 
observing the commutator relation we are left with 

^7?1oc(m) = ■^A([i/int(M - 1, M), Aoc(m)] + [Hinti^i, ^ + 1) , ^loc (A*)] ) ■ (29) 

In sofar as the local energy is a conserved quantity, this equation may be rewritten 
as 

^ilioc(M) = divJ = J(^, ^ + 1) - J(^ - 1, ^) , (30) 
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where we have introduced the discrete version of the continuity equation. Thus, the 
right hand side of ([29|) can be interpreted as a current into as weh out of the central 
subunit /i and, thus, we define the current operator as 

J(/i,M+l) = ■^A[7?i„t(/i,/i+l),i?ioc(Ai)] • (31) 
3.3. Open Quantum Systems 

The time evolution of the above defined system is governed by the celebrated 
Liouville-von-Neumann equation 

Switching from the Hilbert space Ti, to the Liouville space C it is convenient to write 
the last equation as 

= (33) 

defining £ as a so-called Liouvillian—, a super operator acting on operators of the 
Hilbert space, here the density operator. Transforming the density operator to a 
"vector form" we can write the Liouvillian as a tensor of second order (matrix) 
of squared dimension again (in fact a super operator is usually a tensor of higher 
order) . 

The above Liouvillian only describes a coherent time evolution of the closed 
quantum system. The von Neumann entropy of the complete system remains con- 
stant. However, the Liouvillian can be extended to a more general form containing 
also damping processes due to the coupling to an environment. The environment 
itself is not described as a concrete physical system, but only by its action on the 
system. This coupling leads to equilibration and to a maximum entropy state in the 
considered system, thus, a thermodynamical situation. 

The only restriction on the respective super operator £, whether coherent or not, 
is that the respective dynamics has to map a density operator again to a density 
operator. In 1974 Lindblad^® was able to introduce the most general (Markovian) 
form of such an operator. Thus the time evolution equation for an open quantum 
system is given by the quantum master equation in Lindblad form— 

(£s+£E(r,AE)) /) (34) 

with the coherent part £3 and the dissipator £e- Ae denotes the system bath cou- 
pling strength. In case of a single two-level system {n — 2), e.g., controlled by two 
damping channels with rates Wqi and Wio, depending on the bath temperature and 
the bath coupling strength, we obtain for the dissipator 

£e (r, Ae ) /O = M^io (T, Ae ) (2o-_ pa+ - pcT+ a- - a+ ct- p) 
+ Woi(T, Ae)(26-+P(t_ - p<T-<T+ - <T-o-+/o) , 
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in terms of raising and lowering operators (t_|_, (7_. Using this super operator at 
the edges of a spin chain by simply generalizing the Lindblad operators to raising 
and lowering operators times the unity of the rest of the system one gets a local 
bath coupling. Thus, a heat conduction scenario is implemented (see Fig. [T]). This 



































E 


1 








f 








J 








1 


F 


1 










I 












Ae 








Ae 










T + AT/2 














T- 


AT/2 


Reservoir 1 




System 






Reservoir 2 



Fig. 1. Quantum mechanical iicat conduction model system. Several weakly coupled spins (A) 
coupled to two heat reservoirs (Ag) of different temperature modeled by a quantum master equa- 
tion. 



procedure has been criticized—, but for weakly coupled subunits it is a valid method, 
nevertheless. A slightly different approach has been used in the literatur e 
by diagonalizing the system first, transforming the operators into this eigenbasis. 
However, both approaches lead, at least qualitatively, to the same results. 

3.4. A Heat Conduction Model 



0.2 
^ 0.15 
0.1 



12 3 4 

system number n 



0.2 
<^ 0.15 
0.1 



12 3 4 
system number /i 



(a) Heisenberg chain (b) Forster chain 

Fig. 2. Temperature profile for the model in Fig. [T] (four spins, bars in the boxes mark the 
appropriate temperature of the bath), (a) Heisenberg model (b) Forster model. 



To describe a proper nonequilibrium scenario we need several ingredients: a 
spatially structured system and a local environmental coupling. The complete Li- 
ouvillian of a heat conduction model system, as depicted in Fig. [U thus reads 

£ = £s + £i(ri,Ai) + £2(T2,A2) . (35) 



I 



I 
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This operator describes a spatially structured model system (the first term defines 
the coherent part) coupled to two heat baths of different temperature at either end, 
given by the dissipators £i and £2- 

In case of a Heisenberg spin chain, e.g., we find a normal diffusive behavior, i.e. 
a finite temperature gradient (cf. Fig. [2lja)) within the system and Fourier's Law to 
be fulfilled (cf. Fig. [3|). On the other hand, a Forster coupled chain shows ballistic 
behavior, i.e. a flat temperature profile and thus a violation of Fourier's Law (see 
Fig. [2Kb)). See Refs.— , e.g., for a detailed discussion of the transport behavior in 
open spin chains. 



■a 




0.0001 0.0002 0.0003 

AT(2,3)[A£;] 

(a) Current 




0.1 0.2 0.3 0.4 

r[A£] 

(b) Conductivity 



Fig. 3. (a) Validity of Fourier's Law for Heisenberg model: Current over temperature difference 
between the two central systems (2, 3). The temperature differences of the heat baths are taken at 
mean temperatures T = 0.00015,0.2015,0.4515, respectively, (b) Dependence of the conductivity 
K on the mean temperature. 



4. Perturbation Theory in Liouville Space 
4.1. Liouville Space Algebra 

Following Schack— et al., the set of linear operators acting on states of an n- 
dimensional Hilbert space Ti. constitute a d = n^-dimensional complex vector space 
- the Liouville space C. Sorting all entries of an operator A in Ti. into an n^- 
dimensional vector we can introduce "ket" and "bra" vectors in C 

A^\A) , A^ ^ {A\ . (36) 

As a complex vector space one then defines an inner product in C given by the 
trace-norm of operators in Ti. 

{A\B) = Tt{A1b} . (37) 

Operators acting on these "states" (operators in Ti) in the Liouville space are called 
super operators. A super operator 21 = |y4)(i?| acts on an arbitrary state \X) of C 
according to 

a|X) = \A){I3\X) = Tr{B^X}i . (38) 
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A complete set of operators in H constitutes a basis for the Liouville space 

of the respective system, e.g. the set of Pauh operators {l,ax,(Jy,a-z} for a spin. In 
general this set of operators is complete but not orthonormal, i.e. the super operator 

6 = (39) 

3 

is not the unit operator i in Liouville space. According to the definition of the 
transformation operator in the basis {|^j)} 

(40) 

one can rewrite (|39p as 

= ^|i,)(i,H^J^jt. (41) 

j 

Thus, the determinant of & reads 

det © = det (2J 2J^) det ^ det «»t . (42) 

All basis operators are linearly independent, det^& 7^ 0, and thus dot© ^ 0, too. 
Therefore the operator 03 is invertible and one finds a dual basis \A^) according to 

(43) 

^|i,)(i-'>i. (44) 
i 

detailed introduction to super operators. 
4.2. Unperturbed System 

The model depicted in Fig.[T]is described by the Liouvillian psp. As the unperturbed 
system we choose both reservoirs to be at the same temperature Ti ^ T2 — Te with 
equal coupling strengths Ai = A2 = Ae, defining £0 = £3 + 2,i{Te, Xe)+ ^2{Te, ^e)- 
Thus, the stationary state of the unperturbed Liouville-von-Neumann equation is 
a thermal equilibrium state po with temperature Te: This state should support 
neither a heat current nor temperature gradients. 

The eigenvalues and eigenvectors of the unperturbed system are given by the 
eigenequation £o\pj) = lj\pj) {j — 0, ■ ■ ■ , — !)■ The above introduced thermal 
equilibrium state po is an eigenvector of the system with eigenvalue zero, £o|/5o) = 
(Zq — 0). Since po is the unique stationary state all other eigenvalues have a negative 
real part. Furthermore, all eigenvectors except pQ are trace free^. 



with the property 

3 

We refer to Refs.— i22i22ii2i for a more 
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4.3. Local Equilibrium State 

The system is perturbed now by applying a small temperature gradient AT, defined 
by the super operator 

£'(AT, t) = £i (r + ^/ W) + £2 (t - ^/W) , (45) 
with the "switch-on function" 



fit) = e(-t)e* + e(t) (46) 

where Q{t) denotes the step function. Thus, we start at time t = —00 and allow 
the perturbation to increase exponentially till t = 0, setting AT constant for all 
times t > 0. Thereafter the system is subject to this small constant external tem- 
perature gradient. Finally, we are interested in the properties of the stationary local 
equilibrium state of the system reached in the limit t —>■ 00. 

The time evolution of the whole system under the influence of the perturbation 
is given by the Liouvillc-von- Neumann equation 

||/5) = (£o+£'(Ar,0)|/5). (47) 

Starting from a thermal equilibrium state po in the past, we assume the time de- 
pendent state of the whole system to be 

\m) = \po) + • (48) 

Introducing this into (|47p . suppressing terms of higher order in the perturbation 
and observing that dpo/dt = 2o\po) = 0, one finds 

^ |Ap(t)) - £0 I A/3(t)) = £'(Ar, t) \po) . (49) 



Using the operator transformation in Liouville space as introduced by Kub( 



,15 



l[e~^«^\Ap{t))]) = ||Ap(i)) -£o|Ap(0) (50) 



one finds 

^ [e-^«*| Ap(i))] = e-^«*£'(Ar, t) \po) . (51) 
Formally integrating this differential equation, we get 

\Ap{t)) = /* dt' e^°(*-*') £'(Ar, t') \po) . (52) 
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4.4. Standard Kubo-Formula 

For the special system described by the LiouvilHan 

^ = {£o + 2')p=-^[Hs+m,p] , (53) 

one could derive the standard Kubo- formula from (|52p . where the perturbation 
is considered a potential part H^. in the system Hamiltonian Hs- Here, we allow 
for explicitly time-dependent perturbations H'^ = F{t)H' with an explicitly time- 
independent operator H' and an external field F{t), containing switch on functions 
like ((46|) as well as some oscillating terms with frequency lo. Note that also the 
unperturbed system does no longer contain any thermal reservoir. 

The interesting quantity in context of nonequilibrium states is the heat current. 
Due to the definition of the current operator in Sect. 13.2] and the fact that po does 
not support any currents one gets 

J = Tr{ Jp(t)} = Tr{J/5o} + Tr{JAp(t)} = Tr{JAp(i)} . (54) 

Finally, considering the special system ([SS]) one arrives at the standard Kubo 
formula?^ 

/•oo rP 

J^j dt'F{t'~t) d(3'Tr{J{~t' ~ihl3')Jpa} , (55) 
Jo Jo 

where the first current operator refers to the Heisenberg picture (cf. (|23p ). Basi- 
cally, this is a current-current autocorrelation function for the linear response of a 
quantum system to an external perturbation. 

The above method to investigate the perturbation of a quantum system is re- 
stricted to perturbations which are defined as a term in the Hamiltonian of the 
system. This is not the case for thermal perturbations. These are defined in the 
Liouville space rather than in the Hilbert space. Furthermore, it remains unclear 
how the system reaches the thermal equilibrium state pa without any heat bath 
present. 



4.5. Kubo-Formula in Liouville Space 

We return to the perturbation superoperator (j45|) . This time-dependent superoper- 
ator— could be written as 

2'{AT,t) = ^^fm (56) 

with a time independent perturbation operator € and a time dependent prefactor. 
Thus, ^ reduces to 

\m)) = / dt' e^''(*-*')/(i') e \po) . (57) 
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Introducing the unit operator (|44p in the eigenbasis of the unperturbed system 
yields 

|A/5(i))-^ r dt'e^°(*-*')^|p,)(P^I/(0«|Po). (58) 

With e^"(*-*')|pj) = e'^(*-*')|pj) we find 
A.T\ 

\mt)) = ^r^Y.(P'\^\P^)\P^) / e'^-^*-*')/(i')dt' • (59) 

By integrating over t' this reduces to 

|Ap(t)) - (l^ + ^)(P^I^IPo) \P.) ■ (60) 

This perturbative term wiU include aU currents and local temperature gradients of 
the system under the given perturbation. 

Since we are interested in a local equilibrium state - a stationary state with 
a constant current and temperature profile, which will be reached after a certain 
relaxation time - we consider IjGOp in the limit oi t ^ oo finding 

|Ap) = lim \m)) = .^Y.^mPA 1^^.) . (61) 

This is the first-order change of the density operator introduced by the perturbation. 
Let us call this equation the Kubo-formula in Liouville spaced. Note the similarity 
with the Hilbert space perturbation theory: the change of the density operator 
due to the perturbation depends on both the matrix element of the perturbation 
operator and all eigenstates of the unperturbed system. 



4.6. Heat Transport Coefficient 

Now we are prepared to account in such Kubo-Liouville scenarios for the local 
temperature profile and the expectation value of the current. The stationary density 
operator of the system is given by \p) = \pa) + |A/5). Since we know that jpo) does 
not contribute to any local temperature difference or heat current, the respective 
expectation value is determined by |A/5) only. The local energy difference is given 
by the operator A_ff(/i,/i + 1) = H\oc{fJ.) — H\oc(p + 1), and therefore we find for 
the local internal temperature gradient 

ST{,, A. + 1) = E Tr{A7?,e(M, M + l)p,} (62) 

and the local current within the system 

J(M, M + 1) = E Tr{ J(M, M + Dp,} . (63) 
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The current as well as the local temperature gradient are thus found to depend 
linearly on the global temperature difference of the bath systems. Under stationary 
conditions the current must be independent of /i, J(/i, /i + 1) = J, so that (|63p can 
be rewritten as 

J = -k'AT , (64) 
defining the so-called global conductivity as 

,,^h^Y.imMTr{Ji,,, + l)p,}. (65) 

j=i J 

Eigenstates and eigenvalues entering this global conductivity k' only depend on the 
mean temperature of the unperturbed system, but not on AT. Based on this K as a 
global property of the system including its contact properties, we call (j64l) external 
Fourier's Law, for which the global conductivity (|65p defines the overall resistance 
of our given open quantum system. 

Furthermore, combining (|62p and (|63p . we define a local conductivity as 

implying K{^,fi+ 1) to be independent of the external difference AT. 

This new approach does not have the problem of introducing a potential term 
into the Hamiltonian of the system, like in standard Kubo formulas for heat con- 
duction. The bath systems, modeled by a quantum master equation, directly define 
the perturbation in Liouville space. Like in standard perturbation theory in Hilbert 
space the first order correction to the stationary state of the system is expressed 
in terms of transition matrix elements of the perturbation operator and the eigen- 
states and eigenvalues of the unperturbed system. Only the non-orthogonality of 
the eigensystem of the unperturbed Liouvillian calls for a more careful treatment; 
formally the equations are very similar. 



5. Hilbert Space Average Method for Heat Conduction 

As noted already by Einstein^i^-^ in his study of Brownian motion, there is a close 
connection between the bath assisted transport as considered in the last Section 
and the decay from a nonequilibrium initial state into the global equilibrium within 
a closed system. An important role should be played by the transport coefficients as 
material constants. However, in small closed quantum systems we do not find a regu- 
lar decay behavior, but rather coherent oscillations. Unfortunately, there are neither 
theoretical nor numerical investigations of larger systems feasible at the moment. 
Thus, we turn to yet another approach to heat transport within quantum mechan- 
ics based on the Hilbert Space Average Method (HAMl^^^^^^LS^) . This method 
allows for a prediction of the Schrodinger evolution of certain "coarse grained" ob- 
servables, such as, e.g. the occupation probability of a whole energy band, and is 
thus well suited for the investigation of the decay in mesoscopic quantum systems. 
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5.1. Mesoscopic Model System 

The class of systems we are going to analyze next is depicted in Fig. lU again charac- 
terized by the Hamihonian (|25|) . Instead of using two level systems the N identical 



i 




^J. = l n = 2 ^ = N 



Fig. 4. Heat conduction model: A'^ coupled subunits with ground level and band of n equally 
distributed levels. Black dots refer to used initial states. 



subunits now feature a non-degenerate ground state and a band of n excited states 
each, equally distributed over some band width 5e in such a way that the band 
width is small compared to the local energy gap IS.E {5e ^ AS, 5e in units of 
Ai?). These subunits are coupled by an energy exchanging next neighbor interac- 
tion i/int(Ai,M+ 1); chosen to be a (normalized) random Hermitian matrix allowing 
for any possible transition such as to avoid any bias. Our results will turn out to 
be independent of the exact form of the matrix. We choose the next neighbor cou- 
pling to be weak compared to the local gap (A ^ Ai?, A in units of Ai?). This 
way the full energy is approximately given by the sum of the local energies and 
these are approximately given by ('0(t)|i?ioc(/i)|'0(^)): where is the probability 
to find the /i-th subsystem in its excited state. This model is primarily meant to 
demonstrate how energy transport might emerge from Schrodinger dynamics - a 
statistical diffusive behavior from a time reversible microscopic theory. 

Restricting ourselves to the single excitation subspace, i.e. initial states with 
one system somewhere in the upper band, all others in their ground state, we can 
derive a reduced Hamiltonian model. The respective Hamiltonian matrix is shown in 
Fig. [51 here for three subsystems only. The dark gray blocks refer to the interaction 
Vij between adjacent subunits, whereas the light gray ones to the local Hamiltonian. 
We skip all diagonal couplings, transitions within the band, for simplicity. 

We introduce projection operators P^. These operators are matrices with the 
unit operator of dimension n at position (/z, /i) and zeros elsewhere. Thus, the 
operator projects out the part of the state describing the /ith system being 
in the upper band. The probability of finding subunit /i excited is thus given by 
P^ = Tr{|7^)(V'|P4. 

6. Hilbert Space Average Method 

The Hilbert space average method (HAM) is in essence a technique to guess the 
value of some quantities defined as a function of \ip) if itself is not known in full 
detail, only some features of it. Here the method is used to guess the expectation 
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Fig. 5. Reduced Hamiltonian (local part light gray, interaction dark gray) of the one excitation 
subspace in order to investigate heat conduction in the model displayed in Fig.|4]{JV = 3 subunits). 

value (t/'IAIV') if the only information about \^p) is the set of expectation values 
{il;\Pf^\tp) = Pf^. Such a statement naturally has to be a guess, since there are in 
general many different that are in accord with the given set of P^, but produce 
different values for {tp\A\tp). Thus, the key question for the reliability of this guess 
is whether the distribution of the {ip\A\tp) produced by the respective set of \'ip)'s is 
broad or whether almost all those yield {tp\A\ijj) that are approximately equal. 
It can be shown that if the spectral width of A is not too large and A is high- 
dimensional almost all individual \^) yield an expectation value close to the mean 
of the distribution of the (V'l^lV') a-nd thus the HAM guess wiU be rehable^. To 
find that mean one has to average with respect to the \tp). This is called a Hilbert 
space average A and will be denoted as 

^ = Tr{ia} with a:=[|^)(VlI{(v,|P,|^>=P,}- (67) 

This expression stands for the average of (i/ijAl?/;) = Tr{yl|'0)(V'|} over all that 
feature (-^iP^j^) = but are uniformly distributed otherwise. Uniform distribu- 
tion implies invariance with respect to all unitary transformations U that leave P^ 
unchanged, i.e., {tp\W Pf^U\ip) = P^. Thus the tj are specified by \(j,Pfj\ = 0. 

Any such U has to leave a invariant, i.e., [t/, d] = 0. Furthermore, a has to obey 
Th:{6iP^} — P^. These conditions^i^^ uniquely determine a as 

« = E — (68) 

i.e., an expansion in terms of the projectors P^ themselves. 

Consider a pure state of the total system at some time t, \ip{t)). In the interaction 
picture the dynamics of the full system is controlled by the interaction V{t), only. 
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The time evolution is generated by the corresponding Dyson series truncated to 
second order in a weak interaction scenario 

mt + r)) « (i - ^Uiir) - ^U2ir)) m)) = ^Wl^W) , (69) 
with the two time evolution operators 

Ui{t) = r dT'V{T') , U2{t) = r dr' r dT"V{T')V{T") (70) 



and the interaction Hamiltonian V(t) in the interaction picture. This allows for the 
computation of the probabilities at time t + r finding the /ith system somewhere 
in its excited band 

P^it + t)= Tr{D{T) \m){m\ D\r) . (71) 

Assume that rather than \tj}{t)) itself only the set of expectation values P^ is known. 
The application of HAM produces a guess for the P^{t + r) based on the P^{t) 
replacing |?/;(i)} (V'(i)l by the above d ([55]) 

P,it + r) « Tr{^(r) a D^r) P,} ^Y.^ Tr{I)(r) P, D^r) P,} . (72) 

Using this yields after lengthy, but rather straight forward calculations 

P^{t + r) - P^{t) « /(r) (P^_i(t) + P^+i(i) - 2P^(t)) , (73) 

with 

fir) - r r Tt{V{t") V{0)} dr" dr' . (74) 
nJo Jo 

The above integrand is essentially the same environmental correlation function that 
appears in the memory kernels of standard projection operator techniques. Those 
correlation functions typically feature some decay time Tc after which they vanish. 
Integrating twice yields functions that increase linear in time after Tc. Hence, for 
T > Tc one simply gets /(r) — ^t, where 7 has to be computed from ([74)) . but 
typically corresponds to a transition rate as obtained from Fermi's Golden Rule, 
i.e., 

Assuming that the decay times of the correlation functions are small one can trans- 
form the iteration scheme (j73p into a set of differential equations 

d P 

^ = -7(^1-^2), (76) 
d P 

^ = -7(2P,-P^_i-P^+i), (77) 
^ = -7{Pn - Pn-i) . (78) 
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The above scheme, however, only applies if the dynamics of the system are 
reasonably well described by a Dyson series truncated at second order also for 
times r larger than Tj,. This and similar arguments yield the following necessary 
conditions for the above described occurrence of diffusive transport 

2a£>1, A2^«1. (79) 

If those criteria are violated no diffusive transport is expected. For a more detailed 
description of HAM and its various implications, see Gemmer— 

6.1. Complete Solution 

To analyze validity and performance of HAM we compare its results with data 
from a direct numerical integration of the Schrodinger equation. This is, of course, 
only possible for systems small enough to allow for the latter. Hereby we restrict 
ourselves to initial states with only one subsystem in the exited band (all others in 
the ground level, black dots in Fig. [4]). Finding an effective Hamiltonian for the one- 
excitation subspace we are able to solve the Schrodinger equation for up to TV = 10 
subsystems, n = 500 levels each. Firstly restricting to = 3, the numerical results 
together with the HAM predictions for an initial state with Pi = 1, P2 = ^3 = 
are shown in Fig. Eta) (A^ = 3, n = 500, Se = 0.05, A = 5 • 10"^). There is a good 
agreement. 




(a) Decay (b) Accuracy 

Fig. 6. (a) Probability to find the excitation in the 11= 1,2,3 system. Comparison of the HAM 
prediction (lines) and the exact Schrodinger solution (dots). {N = 3, n = 500, A = 5 ■ 10~^, 
Se = 0.05) (b) Deviation of HAM from the exact solution (Df least squares): dependence of Df 
on n for N = 3. Inset: dependence of Dj on A'^ for n = 500. 

To investigate the accuracy of the HAM for, e.g., Pi{t) we introduce Df, being 
the time-averaged quadratic deviation of the exact (Schrodinger) result for Pi{t) 
from the HAM prediction 

Dl ^ 1 J\p^^^\t) - P^^'^'^^t jfdt . (80) 

To analyze how big this "typical deviation" is, we have computed for a A^ = 3 
system with the first subunit initially in an arbitrary excited state, for different 
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numbers of states n in the bands. As shown in Fig. [Hb), the deviation scales hke 
1/n with the band size, i.e., vanishes in the hmit of high dimensional subunits. 
This behavior does not come as a complete surprise since it has theoretically been 
conjectured and numerically verified in the context of equilibrium fluctuations of 
non-Markovian system s The inset of Fig. [6ljb) shows that goes down also 
with increasing number of subunits N but then levels off. Altogether HAM appears 
to be applicable even down to moderately sized systems. So far we have restricted 
ourselves to pure states. A drastic further reduction of Di can be expected for mixed 
states (which are typical in the context of thermodynamical phenomena), since pure 
state fluctuations should cancel partially if added together. 

6.2. Energy Transport 

In order to account for the energy diffusion constant in the given model system, 
we have to consider the energy current inside the system. The total internal energy 
in one subunit of the above described type is given by the probability to be in 
the excited band of the subunit times the width of the energy gap, = AEP^. 
Thus the current is defined by the change of the internal energy of two adjacent 
subunits 





fdUl 


dU2\ 




fdPi 


dP2 


I 


{ dt 


dt J 




{ dt 


dt 



(81) 



The change of the probability in time is given by the rate equation (|76p and one 
finds the current 

J = -~kAE{P2 - Pi) - -n{U2 - Ui) . (82) 

Thus, the current is a linear function of the probability gradient respectively energy 
gradient inside the system. This specifies a diffusive behavior since energy spreads 
through the system in a statistical way. Therefore we may extract the energy diffu- 
sion constant from the above equation, identifying the parameter k, from the rate 
equation with the transport coefhcient given in (j75p . This energy diffusion constant 
depends on the coupling strength of the subunits A and on the state density of 
the excited band n/5e. Remarkably this diffusive behavior is not restricted to an 
initially small energy gradient. 

6.3. Heat Transport 

So far we have considered energy diffusion through the system only. The final state 
should approach equipartition of energy over all subunits - a thermal equilibrium 
stated. Close to this equilibrium we expect the system to be in a state where the 
probability distribution of each subunit is approximately canonical (Gibbs state), 
but still with a slightly different temperature for each site. Specializing in those 
"local equilibrium states" and exploiting the HAM results allows for a direct con- 
nection of the local energy current between any two adjacent subunits with their 
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temperature gradient AT ^ T1—T2 and their mean temperature T = (Ti + T'2)/2. 
Since this connection is found to be hnear in the temperature gradient one can 
simply read off the temperature dependent heat conductivity— 



Me [Tt 



AE 

e 



AE 

ne kT 



(83) 



as displayed in Fig. [71 For the present model this is in agreement with k — jc 
according to (|75[) . if one inserts for c the specific heat of one subunit. 




6q 
<1 

CO 

O 



4 6 
kT[AE] 

Fig. 7. Heat conductivity II83I I over temperature for a system with n = 500, 5e = 0.05 and 
A = 5 ■ 10-5. 



7. Conclusion 

Is the question concerning the microscopic foundation of Fourier's Law completely 
answered now? Certainly not in all details, however, there are new building blocks 
and perspectives accessible so that a satisfying answer in the near future should 
come into reach. 

Both, the Hilbert Space Average Method and the perturbation theory in Liou- 
ville space, have opened new perspectives towards the investigation of statistical 
behavior in small quantum systems near the global equilibrium. Several interesting 
open questions remain which hopefully will be answered in the near future. Those 
questions mainly address the relationship between the presented new methods, the 
Hilbert Space Average Method and the Kubo formula in Liouville space as well as 
the connection to some other important theoretical concepts like the standard Kubo 
formula, the quasiparticle approach, quantum chaos etc. Work is in progress to ap- 
ply the new methods to other concrete models and to discuss differences between 
the theoretical predictions. However, until now the application of the Hilbert Space 
Average Method on models like spin chains remain challenging. 

In face of the rapid miniaturization of many technical devices like e.g. computer 
chips, and the trend to utilize quantum mechanical systems, a deeper understanding 
of transport and relaxation processes and the thermodynamic properties of small 
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systems will become necessary in the near future. At least, it will be indispens- 
able to know to what extent the thermodynamic concepts can be extended to the 
nanoscale and which limits should be taken into account. Besides the foundation of 
relaxation and transport on quantum mechanics also new aspects beyond a proper 
thermodynamical description are interesting. For example, there are some materials 
featuring extremely high thermal conductivities, due to magnetic excitations. Here 
we found that besides the emergence of equilibrium properties in small quantum 
systems, nonequilibrium thermodynamic aspects are already present far below the 
thermodynamic limit. Hopefully, the ongoing research in this field will help to shed 
new light on the fundamental aspects behind those interesting and long since known 
topics of relaxation and transport. 
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